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We study theoretically inelastic spectrum of coherent backscattering of laser light 
by two atoms. For an intense laser field, there are frequency domains of not only 
constructive but also destructive (self-)interference of the inelastic photons. We 



interpret the emergent spectral features using the dressed states and considering 

2 ■ 

\ coherent backscattering as a kind of the pump-probe experiment. 
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I. INTRODUCTION 
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alization of Bose-Einstein condensate |2j made it also possible an exploration of various 



mesoscopic transport and localization phenomena 



usi 



using cold atomic gases. In particu- 



lar, possessing large (compared to their geometric size) and easily manipulated scattering 



cross-section, atoms turned out to be suitable for studying light transport in the weak and, 

3 : 

prospectively, strong localization regimes. Recently, coherent backscattering (CBS) of res- 



onant laser light (an analog of weak localization o 



c 



ectrons in disordered conductors) has 



been observed with cold, trapped atomic clouds "\ Since then, CBS of light by cold 
atoms has become an area of intense theoretical and experimental research (for a recent 
review, see jf]]). 

Coherent backscattering is an enhancement of the average intensity of light reflected off 
a dilute, disordered medium in the backscattering direction. The underlying physical reason 
for the emergence of CBS is the constructive interference between the counterpropagat- 
ing (labelled "direct" and "reversed") multiple scattering amplitudes. When a scattering 
medium consists of individual atoms, several mechanisms affecting phase coherence between 
the interfering amplitudes should be considered. These are (i) Raman scattering on degen- 
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erate atomic transitions; (ii) inelastic scattering; (iii) mechanical motion of atoms. 

As regards the atomic degeneracy mediating Raman processes accompanied by photon 
polarization flips, its dephasing role is nowadays very well understood, with quantitative 
accordance between theory 

ABB 

and experiment p, LH| , for ensembles of Rb atoms. 

The next dephasing mechanism, the inelastic photon scattering induced by atomic satu- 
ration s, has been studied in less detail. Recent experiment with cold Sr atoms demonstrated 
a rapid decrease of CBS enhancement factor a versus saturation at moderate s < 1 [ill ]. 
It was shown, within a scattering theory approach applied to two atoms in the regime of 
weakly nonlinear scattering that a decrease of CBS enhancement factor occurs due 
to the partial distinguishability of the interfering amplitudes. In the general case of many 
atoms, three different amplitudes interfere constructively in the weakly nonlinear regime, so 
that a may exceed the linear barrier 2 1^ . 

In this contribution, we will also be concerned with the impact of saturation on CBS, 
though for arbitrary laser field intensities. But prior to proceed with a presentation of this 
work, let us make a brief note on the effect of thermal motion of atoms at temperatures on 
the order 100 fiK typical for CBS experiments. 

Mechanical motion spoils phase coherence between the direct and reversed amp litudes [l^ 
if spread v of the atomic velocities violates the resonance condition kv <C 27 (l5( , where k is 
the wave number and 27 is the natural linewidth of the excited atomic state. In the regime 
of weak laser intensities, this inequality is usually satisfied, and the picture of motionless 
atoms works very well so long as the CBS intensity is concerned Q|. However, already in the 
elastic scattering regime, the photon recoil and Doppler effects do modify the CBS spectrum 
[rI . It is even more so in the inelastic scattering regime, when the atoms from the cloud are 
rapidly accelerated out of resonance by a powerful laser field. Nonetheless, we will ignore 
this acceleration and assume that the atoms are fixed in space. Thus, we focus here on 
the laser field coupling exclusively to the atomic internal degrees of freedom, in order to 
highlight the fundamental interference effect under the influence of the nonlinear scattering. 
Explanation of this influence for atoms at rest is basic to its understanding for atoms in 
motion. 

More specifically, we will study spectrum of CBS by two atoms in the helicity preser ving 
polarization channel. This topic is above all motivated by our previous work 17|, Il8j | . 



where we established existense of the residual CBS contrast in the deep saturation regime, 
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due to the constructive (self-)interference of inelastically scattered photons. However, this 
constructive interference is a net effect of all inelastic photons. The question that naturally 
arises, of what the character of interference is at a given frequency, can only be answered after 
looking at the CBS spectrum. In this work we answer this question for the particular case of 
exact resonance. We demonstrate that, in the saturation regime, there are frequency domains 
where the interferential contribution exhibits not only constructive but also destructive 
interference, and employ the pump-probe analysis and the dressed states representation to 
identify the scattering processes that are responsible for the emergent spectral features. Our 
results agree with those derived within the Langevin equation approach jlfll ]. 

The paper is organized as follows: We start with a brief presentation of our model and 
the master equation approach that we are using. In Sect. IIHI we present results for the 
stationary CBS intensity and enhancement factor, and thereafter for CBS spectrum. In the 
last Section, we conclude our work. 



II. MASTER EQUATION APPROACH TO CBS OF LIGHT BY TWO ATOMS 

A. Model and the main quantity of interest 

Details of our approach are given in Ref. Q|. Here, we will only present its brief outline. 
We consider a model quantum system consisting of 2 identical, motionless atoms located at 
positions ri and r 2 , with the distance r 12 = |i"i — r 2 | being much greater than the optical 
wavelength. The atoms are embedded in the electromagnetic bath of quantized harmonic 
oscillators and subjected to an external laser field of arbitrary intensity [See Fig. da)]. 
Coupling to the bath gives rise to the spontaneous emission from the excited state and to 
the far-field dipole-dipole interaction responsible for exchange of photons, whereas coupling 
to the laser field gives rise to the Rabi oscillations of populations and coherences in the laser- 
driven transitions of both atoms. Although this approach can, of course, be formulated 
for atoms with arbitrary internal structure, we choose the ground states of the atoms to 
be nondegenerate, while the excited state 3-fold degenerate [see Fig. Hfb)]. An important 
parameter describing the effect of a laser field on atoms is the so-called saturation parameter 
s = f2 2 /2(7 2 + 5 2 ), where Q is the Rabi frequency and 5 = oj^ — ujq is the detuning of the laser 
field with respect to atomic resonance. As already mentioned, here we will be interested in 
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FIG. 1: Model of CBS with two atoms, (a) atoms (black dots) are driven by laser light with right 
circular polarization, while CBS is observed in the helicity preserving channel, that is, with flipped 
polarization. Photons in this channel appear as a result of double scattering, g is the strength of the 
far-field dipole-dipole coupling responsible for exchange of photons; (b) internal atomic structure 
corresponding to a J g = — > J e = 1 dipole transition, ojq is the transition frequency, 27 is the 
radiative linewidth, Vl is the Rabi frequency. Sublevels |1) and |3) have magnetic quantum number 
m = 0. Sublevels |2) and |4) correspond to m = —1 and m = 1, respectively. Thick solid arrow 
shows laser field driving |1) <-> |4) transition, while dashed arrow shows CBS field originating from 
|1) «-> 1 2) transition. 



how the spectrum of CBS behaves as a function of s at 5 = 0. In this case, all results can 
be deduced in an analytic form. 

Raman processes which can strongly affect CBS do not take place on the J g = — > J e = 1 
transition under consideration. Furthermore, incoherent single scattering contribution can 
be filtered out by looking at CBS, e.g., in the helicity preserving (h || h) polarization 
channel. Precisely this channel was probed in a recent experiment with cold Sr atoms 11 1. 
We consider the particular case of the laser light with the right circular polarization, that 
is, El = e+i, in the helicity basis notation. Hence, CBS with preserved helicity corresponds 
to the nipped polarization e = e_i as shown on Fig. ^ 

Spectrum of CBS to be addressed in this paper is derived from the average value of the 
first-order field temporal correlation function [201 ]: 



G«(r,t;r,tf) = <Tr {p[e ■ E<->(r,t)][e* • E«(r,tf)]}> 



conf 



(1) 



where p is the initial density operator of the atom-field system, / + ) (r, t) is the nega- 
tive/positive frequency component of the electric field operator of the scattered field, and 
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(. . .) conf denotes configuration averaging. The components of the scattered field are the 
retarded fields radiated by the atomic dipoles, 

2 

4neoc 2 r 



E(+) ( r ' t ) = l£k:E D «^) e ~ ik - rQ ' ( 2 ) 



a=l 

where eo is the permittivity of the vacuum, D Q = — e_i<7" 2 + e cr" 3 — e + icr" 4 , with crj^ = 
\k) a (l\ a , is the dipole lowering operator, and t a = t — |r — r a \/c. In writing Eq. (J2J), we have 
assumed that r 12 <C r, that is, the field is detected in the radiation zone at the distance 
much larger than the interatomic distance. In the following, we will for brevity omit the 
r-dependent prefactor of Eq. (j2J) and, consistently, of the temporal correlation functions. 
Inserting Eq. (j2J) into Eq. we obtain, in the steady state limit t — > oo, 

2 

Gg\r) = <<<^)>sse^> conf , (3) 

a,P=l 

where "ss" stands for steady state, r = t' — t > 0, the inner angular brackets indicate 
quantum mechanical expectation value [see Eq. (0)], and r^g = r a — Yp. 
Spectrum can be obtained via Laplace transform of (JHJ) 21]: 

S(v) = ±hmRe{Gg\z)}, (4) 

where Gil\z) = / °° drexp(— zT)Gi^ (t), z = T — iv with T > and u = uj — uj l . Note that 
the spectrum is defined with respect to the laser frequency which means that the atomic 
correlation functions must be evaluated in the frame rotating at u>l- 

Let us conclude this subsection with a remark on the configuration averaging procedure. 
This procedure is necessary because the two-atom correlation functions may sensitively 
depend on the interatomic distance and orientation of the vector r 12 with respect to k^, 
exhibiting rapid oscillations around the backscattering direction. These oscillations have the 
same nature as a speckle pattern scattered off a disordered medium. After many realizations 
of the disorder, all peaks except the one, corresponding to CBS, disappear. A simple and 
sufficient way to mimic disorder in a two-atom system is to assume an isotropic distribution 
of the radius- vector connecting the atoms and a uniform distribution of interatomic distances 
around the average distance I equal to the scattering mean free path. 
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B. Master equation 

To find the atomic correlation functions appearing in the right hand side of Eq. © we 

nn □ 

have adapted [TO, Us] a theoretical approach initiated by Lehmberg in 1970 [22j. Within 
this approach, dynamics of the dipole operators' expectation values as well as dipole-dipole 
correlators is governed by the master equation 

2 2 

(Q) = J2(C a Q)+ ^Q), (5) 

a=l a^/3=l 

where the Liouvillians C a and generate the time evolution of an arbitrary atomic oper- 
ator Q, for independent and interacting atoms, respectively. Explicitly, 



C a Q = -i*[Dt .D ,Q]-i[fi a (Dt . £i ) + ^(D Q . £ *),Q] 

+7 (Dt • [Q, DJ + [Dt , Q] . D a ) , (6) 
£ a/3 Q = Dt ■ ¥ (g, n) • [Q, D^] + [D* , Q] • TT *(g, n) • D a , (7) 

where f2 Q = fie 4 ^' 1 " 01 . The radiative dipole-dipole interaction due to exchange of photons 
between the atoms is described by the tensor T (g, n) = <yg A, with A = 1 — fin being 
the projector on the transverse plane defined by the unit vector n along the connecting 
line between the atoms a and j3. This interaction has a certain strength depending on the 
distance between the atoms, via 

with ko = u>o/c, and on the life time of the excited atomic levels, through 7. The coupling 
constant \g\ <ti 1 is small in the far-field (k r a/ 3 ^> 1), where near-field interaction terms of 
order (/cor^) -2 and (/co r «,a) _3 can be neglected. 

Of course, an arbitrary operator Q inserted into Eq. (j5J) does not result in a closed 
differential equation. Our system consisting of two 4-level atoms leads to 255 = 4 2 ■ 4 2 — 1 
linear coupled equations of motion for the one-time averages. We solve them perturbatively 
up to g 2 , to account for the lowest order (double-)scattering process giving rise to a nontrivial 
interferential contribution. To help the reader keeping this in mind we will supply symbols 
denoting double scattering intensities and spectra with the subscript "2" . 
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Note that Eq. (jSJ) describes evolution of the expectation values (one-time correlation 
functions), whereas G[Q (t) is the two-time correlation function. By virtue of the quantum 
regression theorem |2l| , the latter satisfy Eq. (JHJ) also, but their initial conditions are ex- 
tracted from the stationary solution of In particular, the double scattering counterpart 
of Gss^O) is nothing but the stationary average backscattered light intensity which will be 
referred to as It^. There is an obvious relation between /* ot and S<>.(v): 



oo 



P 2 ot = / duS 2 (u). (9) 

J — oo 

The expression for 1^ can be obtained independently from Q. We will use this independent 
derivation as an implicit verification of our results for CBS spectra. 

The total CBS intensity at the backscattering direction can be decomposed in the sum 
of two terms 

J*°t = L tot + d°\ (10) 
where C l 2 ot = C\ ot (e = 0) (i.e., k = -k L ), and 

C?\6) = 2Re((a 2 1 1 ^ 2 )I 2 s ] e ^) conf , (11) 
LT = (K 1 2 )I 2 s ] + K 2 2 ) [2l )ss)conf, (12) 

are the so called crossed and ladder terms, respectively. Using these terms we can derive a 
standard measure of phase coherence between the counterpropagating amplitudes in CBS - 
the enhancement factor 



a 



tot 
2 



« = i + t4t> ( 13 ) 

which for perfect two- wave interference is equal to 2. 

Generally, the total backscattered light intensity has the elastic and inelastic counterparts, 



/* ot = If + J^ el . (14) 

The elastic counterpart is given by the product of the expectation values of the atomic 
dipoles, 

2 

If = E («)ss(^2)sse- jk -^) CO nf, (15) 
a,P=X 
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wherefrom for a = f3 we obtain the elastic ladder term Lf and for a ^ (3 the elastic crossed 
term Cf . Given and If we can find the fluctuating part of the dipole correlation 
functions defining I™ el . 



III. RESULTS 



A. CBS intensity and enhancement factor 

The interferential contribution and the incoherent sum, Eqs. (II lj) and (fT2"|) . yield results 



2Re{(a 2 \a 1 2 2 )I 2 s ]e^} = \g\ 2 \ A +1 , +1 | 2 ^ 

xcos{(k + k L )-r 12 }, (16) 

(^ 2 )i 2 s ] + k 2 2 )I 2 s ] = i^riA +1>+1 | 2 ^. (ir) 

-Ri(s), R,2{s), and -P(s) are polynomial expressions in the on-resonance saturation parameter 
s = fi 2 /2 7 2 , 

R^s) = - (6912s + 3168s 2 
9 

+264s 3 + 20s 4 + s 5 ) , (18a) 

R 2 ( s ) = - (H52s + 528s 2 + 132s 3 + 7s 4 ) , (18b) 
o 

P(s) = (l + s) 2 (12 + s)(32 + 20s + s 2 ), (18c) 

and A +lj+ i = e +1 ■ A ■ e +1 . 
The configuration average of (fl"6|) and ([TTjl leads to the final result 

° 2 W - (4 + s)P(s)\ 15 35 > 1 Jj 

_ 2|^| 2 P 2 (s) 

L2 " 15P(s) ' (20) 

with g = g\r a a=i- The scattering angle 9 = 2arcsin{|k + kkl/2/^} <C 1 with respect to the 
backscattering direction was assumed to be sufficiently small herein. 

The enhancement factor a(s), Eq. (fT3*j) . deduced from Eqs. (fl~9"j) and (j2TJj) reads 

»M = 1+ (i£&- (21) 
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FIG. 2: Enhancement factor in the h II h channel versus saturation s. Decrease of a at small 



s is described by the linear function 2 — s/4 in accordance with 



Inset describes a in the 



n 



deep saturation regime. Enhancement tends to the limit = 23/21 [18( indicating constructive 
self- interference of inelastic photons. 

and a(0) = 2.0 in the weak field limit, as expected. The dependence of a on the saturation 
parameter is shown on Fig. |21 For small s, enhancement linearly decreases as 2 — s/4, in 
full agreement with the diagrammatic theoretical result [l2;] and in qualitative agreement 



with the result of Sr experiment 11 1. When s increases further, a monotonically drops to 
an asymptotic value Hindoo a(s) = = 23/21 [3| which is strictly larger than unity, 
implying a nonvanishing residual CBS contrast in the limit of large injected intensities. We 
will next show that this residual enhancement is due to inelastic photons only. Indeed, we 
obtained the following result for the elastic ladder and crossed terms 

2\~g\ 2 s 



rel 
L 2 



cf 



(22) 



15 (1 + s) 4 

As seen from Eq. f!22j) . the elastic component shows perfect contrast for all s. In particular, 
it is this component that results in enhancement a = 2 for very small s —>■ 0. However, 
in the deep saturation regime, this component decreases as s~ 3 , while the counterparts of 
the total intensity, Eqs. (JTHJ), (f2T)|) . as s _1 . Herefrom follows our conclusion about the origin 
of the residual enhancement in the deep saturation regime. Explicitly, the inelastic crossed 
and ladder terms obtained by elementary substraction of Eq. (|2*2*|) from Eqs. (fT9"j) and (|2*D*|) 
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read 

incl _ 2\g\ 2 20736s 2 + 23424s 3 + 7108s 4 + 601s 5 + 44s 6 + 2s 7 

2 ~ IE 9(l + s) 2 (4 + s)P(s) ' ( ' 

inel 2\g\ 2 2016s 2 + 2244s 3 + 796s 4 + 146s 5 + 7s 6 

Ls " ~L5 3(1 + s)2P(s) • ( } 

It is easy to verify that lim^^ C™ el /L£ el = 2/21 =««,-!. 



B. CBS spectrum 

Double scattering spectrum of CBS has the elastic and inelastic components. The elastic 
spectrum at the backscattering direction reads 

J?(i/) = I*6(v), (25) 

where 5(u) is the Dirac's delta-function, and if = Lf + Cf, with the ladder and crossed 
contributions defined in Eq. (|22|1. 

Inelastic spectra of the normalized ladder and crossed terms, for increasing values of Rabi 
frequencies, are shown on Fig. El Normalization is chosen such that integrals of L^i^/L 1 ^ 1 
and C , 2 nel (^)/L 1 2 nel over v yield unity and C^/L 1 ™ 1 , respectively. In the deep saturation 
regime, the value of the latter integral tends to the asymptotic value of the interference 
contrast of CBS, — 1. 

In describing the CBS spectrum, it is natural to use two parameters Q and 7 defining 
positions and linewidths of spectral peaks rather than a single saturation parameter s. We 
will utilize s only to check consistency of our expressions for spectra with the results for the 
inelastic intensity and the enhancement factor. 

As seen from Fig. El inelastic spectra are symmetric with respect to the laser frequency, 
for all Q. For a small value of the Rabi frequency Q = O.I7 [Fig. Efa)], spectra of both 
the ladder and crossed contributions have single peaks at v = 0. Interferential contribution 
C^iy) is positive, though C™ el {v) < L l 2 eX {v) indicating that interference is not perfect in 
this weakly inelastic regime. We can derive analytical expressions for the curves of Fig. 0(a) 
by leaving the leading-order contribution to inelastic scattering ~ (fi/7) 4 , corresponding to 
two-photon processes, and neglecting the higher-order terms. 

The ladder and crossed terms yield the compact expressions (henceforth, we will omit 
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FIG. 3: Normalized inelastic spectra of the ladder (solid line) and crossed (dashed line) terms at 
exact resonance for different values of the Rabi frequency: (a) f2 = O.I7; (b) f2 = 7; (c) £1 = IO7; 
(d) n = IOO7. 



the common prefactor 2|g| 2 /15): 

1 /ft\ 4 7 3 (2 7 2 + z/ 2 ) 



: inel 1 



V) 



1 fVL 



T 



(26) 



7TV7/ 2( 7 2 + i/ 2 ) 3 ' ~ 2 v " y (7 2 + ^ 2 ) 3 ' 

It is easy to establish that the expressions in Eq. ()2rij) are consistent with the behavior of 
the enhancement factor in the two-photon scattering regime. Integrating L 1 ^ 1 ^), C™ 1 ^) 
over all frequencies, we obtain the following inelastic ladder and crossed terms for small Q: 

4 nd * /_°° dvLT\v) = (") 4 > ^ - ^ ^Q nc V) = I (") 4 • (27) 

Rewriting Eq. (J27)l in terms of s and combining it with the small-s expression for the elastic 
ladder and crossed terms = Cf ~ s, we recover the linear decrease 

s 



1 s + 3s 2 /2 n 

rv = 1 -l - — ~ 2 

+ S + 7s 2 /4 4" 



(28) 



As the Rabi frequency fl increases further on [see Fig. |3fb,c,d)], qualitative differences 
emerge [on Fig. EKc,d)] in the behavior of the ladder and crossed terms. First, the spectra 
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split into several distinct peaks. Second, the crossed term becomes negative for a range of 
frequencies beyond the central peak. This is a manifestation of destructive self-interference 
of inelastically scattered photons. Note that a similar effect of antienhancement was reported 
for linear double scattering from atoms with Zeeman-shifted hyperfme ground levels. 
New spectral features are robust and become well-separated in the asymptotic limit of intense 
driving [for an example at O = IOO7, see Fig. 0(d)]. We will next address the approximate 
analytic expression for CBS spectrum at O ^> 7, derived in the leading order ~ (7/O) 2 . 

In this case, the explicit expressions for the ladder and crossed spectra can be represented 
by using a function of two real variables x\ and x 2 : 

I X\ 

£(xi,x 2 ) = ~ x 2 + x 2 - ( 29 ) 

Let us mention the properties of £(xi,X2) that are important to us: (i) if x\ = Const, 
then function (j2S| represents a Lorenzian with width x\ and a resonance at x 2 = 0; (ii) if 
X2 = Const, then (j2T)j) describes a resonance of the dispersive type at x\ = 0. 
With the help of the function (1211) , the ladder and crossed spectra are given by 



LT\v) ^ (^) 2 Qi?(7^) + ^(37^) + ^^(37^-2fi) + £(3 7 ,z/ + 2fi)] 



±^[£(3 7 /2, u-G) + £(3 T /2, v + O)] 
9 

+ 4^(57/2, v - O) + £(5 7 /2, v + fi)] 



14 \ 

+—[£(37/2, v - n/2) + £(37/2, v + n/2)] \ , (30) 

CT\v) ~ (2) 2 Qi?(27,z/) + ^(37^)-^[£(57/2,z/-fi) + £(57/2,z/ + fi)] 
+ ^[£(3 7 , u-2Q) + £(3 T , 1/ + 2fi)n 

+ (^) 3 + 0/2, 3 7 /2) - £(z/ - 0/2, 3 7 /2)], (31) 

where the two terms of order (7/O) 3 are retained because they define dispersive resonances 
of C 2 ncl (is) at v = ±0/2. As seen from Eqs. $Eft and JSJ) as well as from Fig. 0(d), both the 
ladder and crossed terms have 7 resonances, the resonances of both the ladder and crossed 
terms at v — and of the ladder term at v — ±0 being sums of two Lorenzians with different 
widths and weights > 0. The rest resonances of the ladder term are also Lorenzians with 
positive weights. Two resonances of the crossed term at v = ±0/2 have the dispersive line 
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shape and, therefore, have no net contribution to the integrated intensity. Furthermore, 
among the rest five resonances of the crossed term which all are of the Lorenzian type, two 
at v — ±Q have negative weights. Thus, inelastic photons (self-)interfere destructively at 
v = ±f2, yet the overall effect of all inelastic processes is constructive. Note also that in the 
frequency domains where interference is constructive, it is also perfect, as can be concluded 
from the equality between the respective weights of the ladder and crossed terms. 

By performing the elementary integrations of Eqs. (|30J) and (|31|). we arrive at the inelastic 
ladder and crossed terms 

4"'^m 2 , Cr^UTf, ,32) 



3 \QJ ' 9 \n 

which are consistent with Eqs. (|23jl . (|2*3|) and, hence, with a = = 23/21. 

Let us now address the interpretation of the CBS spectrum in the limit of intense driving. 

C. Interpretation 

One can understand the structure of the CBS spectrum from the analysis of CBS as a 
specific realization of the pump-probe experiment [a similar view is held by the authors of 
Ref. jlflj]. In the usual setting of such an experiment j^, an atomic transition is simultane- 
ously subjected to two monochromatic fields: a variable-intensity, fixed-frequency driving, 
or pump, field, and a weak probe field with tunable frequency. For different frequencies of 
the probe field it can be absorbed or amplified depending on the intensity of the pump field. 
This occurs because the pump field leads to the energy levels' shifts and broadenings, while 
the weak field transmission spectrum probes these new resonances; hence the name of this 
technique. 

In our case of CBS with two atoms, an intense pump acts in the |1) <-> |4) transitions 
of both atoms, causing an AC Stark shift of the energy levels. In this case it is instructive 
to treat the laser mode as a quantum system strongly coupled to the laser-driven atomic 
transition J25]. The eigenstates of the laser-atom interaction Hamiltonian for 5 = are the 
dressed states 

I ± (N)) a = (|1, N + l) a ± e ik - r "|4, N) a )/y/2, (33) 

where N and N + 1 refer to the number of photons in the laser mode, and a labels the atoms. 
Spontaneous transitions from the dressed states manifold {| ± (N)) a } to {| ± (N — l)) a } 
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FIG. 4: Scattering processes (a)-(l) and radiative transitions (m)-(p), depicted by solid wavy 
arrows, contributing to CBS spectra in the regime of intense driving. Horisontal lines indicate 
dressed states. In the processes (a)-(l), a photon with either of the frequencies emitted by one 
atom (u>l — Q (a-d); col (e-h); ojl + ^ (i-1)) undergoes the Rayleigh or Raman scattering on the 
dressed states | ± (N — 1)). Diagrams (m-n) show radiative cascades in which CBS resonances at 
lul ± f2/2 appear. Level |2, N) can be populated as a result of a multiphoton scattering process 
with the participation of one doubly scattered photon (depicted by dashed wavy arrows). 

lead to emission of the fluorescence spectrum with three symmetric peaks centered at the 
frequencies ojl — £1, ul, and uj^ + VL known as the Mollow triplet 26(. The Mollow triplet 
emitted by one atom plays a role of the probe for another atom. 

Figure |U illustrates the processes that contribute to the CBS spectrum in h || h channel. 
The left part Fig. EJa-1) shows possible one-photon elastic Rayleigh and inelastic Raman 
processes which photons of frequencies ul — Q (a-d), ujl (e-h), and ul + Q (i-1) undergo on 
the dressed states | ± (N — 1)). The right part of Fig. EJ^m-n) shows radiative cascades in 
the dressed state basis leading to resonances in the CBS spectrum at uj^ ± Q/2. These lines 
appear as a result of the spontaneous transitions from the state \2,N — 1) (note that the 
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atomic state |2) is not affected by the laser field) to states | ± (N — 2)). This is the notable 



Autler-Townes doublet 25, 



271, El emitted from the transition between pairs of states of 



which only one is a member of the laser-driven transition. 

One observes that in all the processes except (b) and (1), several transitions participate in 
the creation of a CBS photon. Phases between the participating transitions can be opposite 
due to difference, e.g., in the initial, intermediate, or final atomic states, leading to the 
negative signs of the interferential contributions around ujl ± Q and ujl ± Q/2. 



IV. CONCLUSION 



Using the master equation approach we have analytically calculated a spectrum of CBS by 
two identical, motionless atoms for the case of exact resonance between the laser and atomic 
transition frequencies. This enabled us to analyse in more detail the previously established 
I7I effect of the constructive self-interference of inelastically scattered photons. 

One conclusion of this analysis is following. The enhancement factor based on the total 
backscattered light intensity is a poor measure of phase coherence between the counter- 
propagating waves in the saturation regime, because it has the important information on 
the character of interference at a given frequency integrated out. At intense driving, one 
should rather use spectrally resolved measurements with a filter whose passband Tf satisfies 
7 <C Tf <C Q. Then, tuning the filter on individual peaks of the CBS spectra, one would ob- 
serve either perfect enhancement (at u = ul', ojl ± 2f2), or antienhacement (at u — ul ± O), 
else no net interference at all (at uj = uj^ ± Q/2). 

Another conclusion is that spectral features of CBS can be qualitatively understood 
from analysis of CBS as a kind of the pump-probe experiment. In the limit of an intense 
driving, when the spectral line-shape looks rather complicated [see Fig.0], the 'pump-probe' 
interpretation in combination with the dressed states approach allowed us to identify the 
origin of all the resonances of the CBS spectrum. Although a detailed explanation of the 
character of interference between the different processes, in which photons from the probe 
field (scattered by one atom) are scattered on the dressed states (of another atom), is 
beyond the scope of this work, the following remark is in order. At exact resonance, some of 
the processes interfere constructively while some destructively, with the overall effect being 
constructive. But it is possible to vary populations of the dressed states and, consequently, 
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the weights of different processes by changing the laser-atom detuning S. Nothing forbids 
the overall effect of inelastic photons to be destructive. In fact, we showed in the previous 
work Q| that for large detuning it is indeed so. We also established in the same paper that 
the ladder term of double scattering becomes negative in h _L h channel in the saturation 
regime. This result can be interpreted as a mere absorption of the probe field. 

A challenging problem for the future would be to generalize the results of the present 
work for moving atoms. 
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